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ABSTRACT 

A plasma instability theory is presented for the prompt radiation from gamma-ray 
bursts. In the theory, a highly relativistic shell interacts with the interstellar medium 
through the filamentation and the two-stream instabilities to convert bulk kinetic 
energy into electron thermal energy and magnetic field energy. The processes are 
not efficient enough to satisfy the Rankine-Hugoniot conditions, so a shock cannot 
form through this mechanism. Instead, the interstellar medium passes through the 
shell, with the electrons radiating during this passage. Gamma-rays are produced 
by synchrotron self-Compton emission. Prompt optical emission is also produced 
through this mechanism, while prompt radio emission is produced through synchrotron 
emission. The model timescales are consistent with the shortest burst timescales. To 
emit gamma-rays, the shell's bulk Lorentz factor must be ;> 10 3 . For the radiative 
processes to be efficient, the interstellar medium density must satisfy a lower limit that 
is a function of the bulk Lorentz factor. Because the limits operate as selection effects, 
bursts that violate them constitute new classes. In particular, a class of optical and 
ultraviolet bursts with no gamma-ray emission should exist. The lower limit on the 
density of the interstellar medium is consistent with the requirements of the Compton 
attenuation theory, providing an explanation for why all burst spectra appear to be 
attenuated. Several tests of the theory are discussed, as are the next theoretical 
investigations that should be conducted. 



Subject headings: Gamma rays: bursts 



1. Overview 



Current observations of gamma-ray bursts place a number of strong constraints on gamma-ray 
burst physics. The measured red shift of z = 0.835, 0.966, and 1.61 for lines in the optical 
emission of gamma-ray bursts GRB 970508 (Metzger et al. 1997a,b), GRB 980703 (Djorgovski 
et al. 1998a, b), and GRB 990123 (Kelson et al. 1999), respectively, and the indirect redshifts 
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of z = 3.5, 5, and 1.096 for GRB 971214 (Kulkarni, S. R., et al. 1998), GRB 980329 (Fruchter 
1999), and GRB 980613 (Djorgovski et al. 1999), respectively, show that gamma-ray bursts are at 
extraordinarily-high redshifts. The power-law gamma-ray spectrum above 511 keV and the rapid 
rise times of the gamma-ray light-curves force one to consider theories with highly-relativistic bulk 
motion in order to avoid thermalization of the gamma-rays through photon-photon pair creation 
(Schmidt 1978; Baring, & Harding 1996). The gamma-ray spectrum can be characterized by the 
E-peak energy (E p ), which is the photon energy of the maximum of the vF v spectrum. Because 
the distribution of E p values is narrow, with an average value of E p ~ 250 keV (Mallozzi et al. 
1996; Brainerd et al. 1999), one suspects that the characteristics of the spectrum are independent 
of the bulk Lorentz factor, which should vary greatly from burst to burst. In some bursts, the 
shape of the gamma-ray spectrum is inconsistent with optically thin synchrotron emission (Preece 
et al. 1998a,b). The shape of the burst spectrum is consistent with being Compton attenuated 
by a high density 10 5 cm -3 ) interstellar medium (Brainerd 1994; Brainerd et al. 1996,1998); 
observations show the x-ray excess predicted by this theory (Preece et al. 1996), and the redshifts 
of ~ 1 to ~ 10 derived by fitting the model to burst spectra are consistent with the values 
measured at optical wavelength (Preece & Brainerd 1999). The success of this theory, its presence 
in every gamma-ray burst, implies that a high density interstellar medium is necessary for a burst 
to occur. 

The observations suggest that the sources of gamma-ray bursts are compact — perhaps a 
several mass black hole, perhaps a supermassive black hole — and that these sources eject mass at 
relativistic velocities in one short event. The observed behavior of the gamma-ray burst arises from 
processes that convert the kinetic energy of the ejected material into electromagnetic radiation. 

The general view that has developed concerning the gamma-ray emission of gamma-ray bursts 
is that it is radiated behind a shock that has developed within a shell moving with a bulk Lorentz 
factor of r ~ 10 3 (Meszaros 1998). The radiative mechanism universally cited is synchrotron 
emission from electrons accelerated to an energy close to the equipartition energy (Tavani 1996). 
The shock itself is collisionless, arising either when the shell runs into the interstellar medium, 
or when the fastest portions of the shell overtake slower portions (Rees & Mesaros 1992, 1994; 
Piran & Sari 1998; Mochkovitch & Daigne 1998). This model has several shortcomings. First, the 
gamma-ray burst spectrum is often harder at x-ray energies than is allowed by an optically-thin 
synchrotron emission model. Second, it provides no explanation for the x-ray excess, for the 
narrow E p distribution, or for the apparent absence of bursts with T < 10 3 , as inferred from 
the absence of photon-photon pair creation and thermalization. From the theoretical side, many 
assumptions are made without theoretical development; among these are that collective processes 
exist that mediate the shock, that a strong magnetic field is generated in the shock, and that the 
energy is efficiently transferred from the bulk-motion of the ions to the thermal energy of the 
electrons. 

It is with these difficulties in mind that a new theory for the generation of the prompt 
radiation in a gamma-ray burst is proposed. The theory is a plasma instability theory. In 



-3- 



this theory, the shell of relativistic material that is ejected from the source interacts with the 
interstellar medium through the plasma filamentation instability, which, in the relativistic regime, 
has a higher growth rate than the two-stream instability. Because the mass ejection event must be 
short, of order the burst durations, and because the mass travels 2cT 2 times the burst duration, 
the interaction is of a thin shell with the interstellar medium. As the shell passes through the 
interstellar medium, the interstellar medium collapses into filaments that contain strong magnetic 
fields and high electron thermal energies. It is the interstellar medium that filaments rather than 
the shell, because the interstellar medium has the smaller density as measured in each comoving 
reference frame. The ions in the filaments remain essentially at rest with respect to the observer, 
while the electrons move towards the observer with a high bulk Lorentz factor. The magnetic fields 
generated through filamentation are strong enough to produce gamma-rays through synchrotron 
self-Compton emission. One finds that the theory places lower limits on both F and the interstellar 
medium density through selection effects, and that these lower limits lead to the conditions 
required by the Compton attenuation theory. The theory implies that there exists a class of burst 
that produces intense and prompt optical and ultraviolet emission, but no x-rays and gamma-rays. 

In this article, I give an analytic development of the theory outlined above. This development 
explores the relevant plasma and radiative processes, deriving the selection effects inherent in the 
theory, and ascertaining the aspects of the theory that provide means of observationally testing 
the theory. In §2, the characteristic ratio of shell density to interstellar medium density and 
the characteristic thickness of the shell are discussed. In §3, the growth of the two-stream and 
filamentation plasma instabilities is examined. The saturation of the filaments and the inability of 
the filamentation instability to mediate a shock are examined in §4. The electron thermalization 
and isotropization and the rest frame defined by the electron component are discussed in §5. The 
radiative processes of synchrotron emission and synchrotron self-Compton emission are examined 
in §6, where the characteristic frequencies and emission rates of each are derived. In §7, the 
radiative timescales for synchrotron and synchrotron self-Compton cooling are derived. The theory 
has several natural selection effects that constrain several of the free parameters in the theory. 
These are discussed in §8. The basic theory and the results of this study are summarized in §9. 
This section also contains some suggestions for observational tests of the validity of the theory. 

2. Characteristic Model Parameters 

The model is of a fully-ionized neutral shell of electrons and protons passing through a 
fully- ionized interstellar medium of electrons and protons. The characteristic physical parameters 
that describe this theory are rii sm , the number density of the electrons or of the protons in the 
interstellar material, n' shell , the number density of the electrons or protons in the relativistic shell, 
and r, the Lorentz factor of the shell relative to the interstellar medium. The quantities with 
primes are measured in the shell rest frame, while those without primes are measured in the 
interstellar medium rest frame. An important parameter in the discussion that follows is the 
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parameter 77, defined as 



n 



'"ism /-, \ 



r;; 



shell 



There is no simple connection between 77, T, and the parameters that describe the energetics 
of the system. In particular, to relate 77 to M, the mass per unit ster radian of the relativistic 
shell, and R, the distance traveled before deceleration, requires a good understanding of the 
radiative transfer and plasma physics of the problem. The value of 77 is therefore treated as a free 
parameter, with the limits on the acceptable values of 77 set by a combination of theoretical and 
observational constraints. The former is set in this section, while the latter is set in §8. 

A specific value of 77 that has a physical significance is the value of 77 expected from the 
continuity equation for a relativistic shock. For r> 1, one has n' shell = rii sm T, and 

1 

Vshock = p ■ (2) 

When 77 < 77 s /j OC fc, the interstellar medium density in the shell rest frame is less than the shell 
density, while it is greater than the shell density when 77 > rj s h oc k- 

A lower limit r) m i n on 77 can be derived through a theoretical argument concerning pressure 
equilibrium with in the shell. If one assumes that the interaction with the interstellar medium 
exerts a force only at the front of the shell, then one can relate 77 to the temperature of the shell 
and R, the distance the shell has traveled from the source when T is a factor of 2 below its initial 
value. This defines a maximum value for the shell density, and therefore a minimum density for 
77, because when the deceleration force is spread throughout the shell, a lower pressure is required 
to maintain static equilibrium within the shell. Defining the dimensionless 4-velocity in the radial 
direction as u r , the deceleration of the shell is given by 

M.cdu r . i 

In this equation, the pressure 7/ is the pressure exerted on the shell by the interaction with the 
interstellar medium as measured in the shell rest frame. The factor of T converts this proper 
measure of force into the radial component of the 4-acceleration. The parameter T is the sum of 
the electron and proton temperatures in the shell in units of energy. Changing the derivative on 
the left-hand side of equation (3) into a derivative in r and setting dr w R and du « T gives a 
maximum density of 

Mc 2 T 

n shell ,max ~ 2R^T ' 

The value of R in this equation is bounded by a minimum distance Ro that is set by the case in 
which the interstellar medium is swept up by the shell. The amount of interstellar medium per 
unit ster radian that must be swept up to change T by a factor of 2 is m k, M./V, so 

*o = f-^V) = 3-94 x ID- pc M^nT^T- ^ , (5) 
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where M.27 is the shell mass per unit ster radians in units of 10 27 gm, rii sm is given in units of 
cm -3 , and T3 is the Lorentz factor in units of 10 3 . With these parameters, a shell subtending 4-zr 
ster radians and having M27 = rii sm = T% = 1 will carry 1.129 x 10 52 ergs of energy, which is the 
characteristic value inferred from the observations at gamma-ray energies. Using equation (5) to 
express R in units of Ro in equation (4) gives 

, _ m p c 2 n ism T 2 (Rq\ 3 ( ) 

^shelLmax c rr \ r> I ' \ > 



6T \ R J 

The ratio of rii sm to n' shell is then 



nism 6T /R_Y 

rjmin — , — 2 r2 d ' \' ) 



n 'shell, m ax rn p c 2 T 2 \R J 

The value of rj m i n is strongly dependent on the distance traveled. For deceleration over R = Rq 
with T = m p c 2 r, which is the temperature found when the interstellar medium is swept up 
adiabatically, one finds rjmin ~ T _1 = r] s hockj as expected for a shock. For lower values of T, one 
has rjmin < Vshock unless R > Rq by a sufficiently large value. For T w m e c 2 T, which is the case 
for the theory discussed below, then one has i] < r] s h oc k for R < 6.7 Rq. Lower limits on 77 for 
T = m e c 2 T and several different values of T are given as functions of R/Rq in Figures 1 and 2. 

The value of 77 defines the light crossing time scale across the shell width. The thickness of 
the shell in the shell rest frame is related to the mass of the shell by 

= m P n 'sheii 1 ' ■ (8) 
Expressing R in terms of Rq through equation (5), one finds 

\m p n iS m) (3) (i?) T1, ^ 



= 4.05 x 10 15 cm m\ 7 n ism T| 



where r/_ 3 = 77/IO 3 . The thickness of the shell in the interstellar medium rest frame is I = l'/F. 
Defining the shell's characteristic timescale t s h e u as I = ct shell, one has 



\ R 



1 



1.35 x 10 2 s^f 7 nr^r 3 ^ (^) Vs. (12) 

For r] of order 1/r, one requires that R = IORq at rii sm = 1 cm~ 3 for t s h e ii = Is. Equation (12) is 
plotted in Figures 1 and 2 as functions of 77 with respect to R/Rq for several values of T and t s h e u. 
Because t s h e ii must be less than or equal to the burst timescales, which are often observed to be 
< 1 s, model values for 77 must be 77 <S 1. It is shown in §8 that the theory provides this limit on 77. 
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Fig. 1: Limits on rj. — The variable rj is denned by equation (1) as the ratio of the interstellar 
medium proper density to the shell proper density. The shell deceleration distance R is plotted 
in terms of Ro, which is defined in eq. (5). Limits on the value of rj are plotted as solid curves 
for the values T = 2 x 10 3 , 5 x 10 3 , and 10 , with the lower limits given by equation (7), and the 
upper limits by equation (2). The values of rj and R that produce the shell thickness timescales 
(eq. [12]) of 0.1s, Is, 10s, and 100s, are plotted from bottom to top as a long-dashed, dot-dashed, 
three-dot-dashed, and long-dashed lines, with a line for each value of T, which terminates on the 
proper boundary for that T. Other free parameters were set to rii sm = 1cm -3 , M27 = 1. 



- 7- 




Fig. 2: Limits on r\. — The same as Fig. 1, except now rii sm = 10 5 cm~ 3 . The dotted-lines at the 
bottom of the figure are for a shell thickness timescale of 0.01 s, with the remaining curves the same 
as in Fig. 1. 
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3. Plasma Instabilities 

Within the reference frame of the shell, one initially has a plane-parallel charge-neutral 
density profile through which a neutral and uniform plasma streams. The shell is assumed to have 
no initial magnetic field. The dissipation of the bulk kinetic energy of a relativistic shell to the 
interstellar medium must be through a plasma instability, because the densities of the shell and 
the interstellar medium ensure that the collision mean free path is much longer than the thickness 
of the shell. For a relativistic plasma, the two relevant instabilities are the two-stream instability 
and the electromagnetic filamentation instability (Davidson 1990). For v <C c, the former has a 
larger growth rate that the latter by a factor of c/v, but for F S> 1, the growth rate of the latter 
becomes much larger than the growth rate of the former. 

The two-stream instability has a maximum growth rate in the relativistic regime of 

1 _3 1 I Aire 2 1 

where uj' p e ism is the electron plasma frequency of the interstellar medium in the shell rest frame. 
Taking T = 10 3 and 

Tiism — 1 cm 3 , one finds a growth rate for the electron two-stream instability 
of 72,, e = 30 s~ x in the shell frame. This growth rate corresponds to distances of travel of 10 9 cm 
for the interstellar medium through the shell. 

The filamentation instability acts on counterstreaming plasmas by creating a magnetic pinch. 
Of the two plasmas, the interstellar medium and the shell, that with the smaller number density 
in it's own rest frame is the plasma that filaments. For gamma-ray bursts, this is the interstellar 
medium when 77 <C 1. 

The growth rates of the filamentation instability is easily derived from the dielectric tensor 
for cold electron and ion streams in the relativistic regime and in the absence of a magnetic field. 
The equation must satisfy (Davidson 1990) 



r 3 jC 2 kl Tjio 2 



0, (14) 



where uj p j is the plasma frequency of component j for the density measured in the observer's rest 
frame, Tj is this component's Lorentz factor, and k± is the wave number perpendicular to the 
velocity vector of the streams. 

If one is considering only a charge neutral stream of electrons and ions moving with Lorentz 
factor T through a charge neutral background, and if one defines the background plasma to be the 
shell, then equation (14) gives a frequency of 



f)2 (. ,12 + lJ 2 ^ 

!2 ^ \~p,e,isrn ' ^p^i^ism , 



l + c~V (u' p l shell + u' p 



12 ^ 
shell j 



(15) 
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In this equation, w peism and w' piism are the electron and ion plasma frequencies of the interstellar 
medium, and uo' p e shell and uj' p i shell are the electron and ion plasma frequencies of the shell, all 
measured in the shell rest frame. The negative value of u/ 2 shows that the wave grows. Because 
'p,e,ism i,ism an d ^'p^shell ^ i °''p,i,sheih one fhids that the growth rate of the filamentation of 



iO. 



the electrons is 

The growth rate is independent of wave length for k± > uo' p e s^u/c, and it is oc k± otherwise. The 
length scale of the filament is therefore given by Xf = kj 1 = c/uj P)eyS h e ii- For ions alone, the growth 
rate, which differs from equation (16) only in the numerator, is given by 

, H.vsm 1 I toe 2 i 

Vi = -l > . =5 , =T= ™o\ ■ (17) 



2,./2 

p,e, shell 



This has the same length scale as the electron filamentation instability, but a lower growth rate. 

The filamentation instability growth rate is faster than the two stream instability by the 
factor of T. For r = 10 s , n>ism — 1 cm ^, one finds an electron filamentation growth rate of 
7j e = 6 x 10 4 s _1 , which gives a length scale of 10 6 cm over which the instability occurs. For ions, 



the filamentation instability grows at a rate that is a factor of ^m e /m p slower, so that for the 
parameters given above, 7^ = 1.4 x 10 3 s _1 . The growth rates of all three instabilities are shown 
in Figures 3 and 4 as functions of V. Both filamentation growth rates are therefore higher than the 
two-stream growth rate for electrons using the parameters derived above for gamma-ray bursts. 



4. Filament Saturation 

The filaments grow until the growth rate of the thread equals the magnetic bounce frequency 
of the particle beam producing the instability (Davidson et al. 1972; Lee & Lampe 1973). The 
bounce frequency, which describes the motion of a particle across the filament through the toroidal 
magnetic field, is given by 



= y eB /mcTx j , (18) 

where Xf is the length scale of the filament, and m is the mass of the particles comprising the 
filament. The size of the filaments is set by the lower limit on the wave number for which growth 
occurs, k±c = oj' pe shell . The ratio of the magnetic field energy density to the particle beam energy 
density in the thread is then [35] 

W' B _ m e n ism _ m e ^ ^ 



W'. 



mn',,, m 



One point to note about this equation is that electron and proton components each generate 
magnetic fields of the same strength. For both components, one has W' B = m e c 2 n isrn r]T 2 , which, 
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Fig. 3: Timescales. — Instability, radiative, and isotropization timescales versus Lorentz factor. 
Plasma timescales are given for an ISM plasma density of rii sm = lcm -3 . The lower two solid 
curves are the linear growth timescales of the electron and the ion filamentation instabilities, t'f e 
and tjj. The dotted curve is the timescale for the electron two stream instability. The short dashed 
curve is the timescale for the electron distribution to isotropize, ti so . The synchrotron cooling 
timescale tsync 

is given by the line with two dots and one dash, while the synchrotron self-Compton 
cooling timescale tcomp is given by the dash and dot line. 
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Fig. 4: Timescales. — Instability, radiative, and isotropization timescales versus Lorentz factor. As 
in Fig. 3, but for rii sm = 10 5 cm~ 3 . 
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for r = rf 1 = 10 3 , gives B' = 0.14 G with n ism = 1 cm" 3 , and B' = 45.36 G with n ism = 10 5 cm' 3 . 
This independence is a consequence of the length scale of filamentation being set by the shell 
electron density. 

The thermalization of the ions within the thread can be estimated by examining the equation 
of motion at saturation. A particle's equation of motion perpendicular to the thread for u x <C T is 
(Davidson et al. 1972) 

$ = ( 2 °) 
Solving this equation of motion using a maximum spatial amplitude of Xf, one finds that the 
maximum momentum perpendicular to the filament is u x = XfCObT/c. Replacing Xf with the 
length scale of the filament, and replacing the bounce frequency with the filament growth rate, 
one finds 

u * = iW T - (21) 

V ism 

If u x <C 1, then the approximate energy that goes into thermalizing the particles in the stream is 
approximately given by W th jW[ sm = u x /F, while if u x 3> 1, it is given by W[ h /W( sm = u x /T. As 
a result, 

W' ism \m V m J 

A point to note is that both of the terms on the right hand side are greater than W' B , so that one 
always has W th > W' B . Equations (19) and (22) are plotted in Figure 5 for m = m p and T = 10 3 
and 10 4 . 

From these equations one sees that the filamentation instability cannot mediate a shock. 
If a shock were present, then rj would be given by equation (2). Placing this equation into 
equations (19) and (22) and setting m = m p , one finds that W B /W( sm srf m e /m p T <C 1 and 
Wl/W[ sm = (m e /m p T) l l 2 <C 1. The energy released through the instability is small compared to 
the kinetic energy of the interstellar medium as measured in the shell rest frame. This contradicts 
the Rankine-Hugoniot equations, and so a shock never arises through the filamentation instability. 



5. Electron Thermalization 

Once ion filamentation is complete, electrons will attempt to come into equilibrium through 
the two-stream instability. This instability will drive the electrons to a distribution that is 
uniformly distributed in energy between the rest frame described by the ions in the filaments and 
the ions in the shell. The rest frame of the electron distribution must preserve both the electron 
charge density and the electron current, since the shell plus interstellar medium is charge-neutral. 
These two conditions define a rest frame for the electrons that has a Lorentz factor relative to the 
interstellar medium rest frame of 

r e = , T + 11 . (23) 
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Relative to this rest frame, the shell is moving with 

_ v f + i 

where the double primes are used to denote quantities measured in the electron rest frame. 

In the electron rest frame, the electron density is given by n" = n ism T e + n' shell F", which can 
be written as 

n': = n' 3heU y/l + r? + 2r,T. (25) 

The electron rest frame has nearly the same Lorentz factor as the shell rest frame as long as 
77 <C T" 1 . For 1 S> rj S> the Lorentz factor of the electrons is T e = ^F/27], which is a factor of 
\/y/2rjT smaller than the Lorentz factor for the ions. The value rj = T" 1 is therefore an important 
transition point for the character of the radiation emitted by the shell, because the Lorentz boost 
of the radiation is smaller than the boost associated with the shell when rj 3> T -1 . On the other 
hand, because 77 < 1, F" = F e > T", so the two-stream instability drives the electron distribution 
to the energy defined by F e . 

The electron in the shell rest frame has a Landau radius that is much larger than the filament 
length scale. Using the definitions for Xf, the gyroradius, and equation (19) for the magnetic field 
strength in the shell rest frame, one can write 

— = u e — i] ■ (26) 
Xf m p 

The gyroradius is therefore much larger than the filament width when 

u e »^V (27) 
m e 

For electrons thermalizing to the energy F e , u e « F / y/T+~2r]F, so that equation (27) becomes 

(^n, ifr ? <r- 1 ; 

T» <^ } m \2 , , (28) 

) 2 i^J r] 3 , otherwise. 
One finds that the inequality holds for the upper term in equation (28) whenever 



T > ^Jm p /m e = 42; below we show that F ^ 10 , so if equation (28) is to fail, it will be 
for the lower term. For r)F > 1, the inequality holds as long as 




(29) 

When inequality (29) holds, each electron passes through many filaments in a single orbit, and 
the motion of the most energetic electrons will be as a single fluid, but when the inequality fails, 
the electrons will be confined to the local conditions within each filament, and the motion of each 
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electron will be determined by these local conditions. For the remainder of the paper, we assume 
that equation (29) holds. 

For the electrons to flow through the magnetic fields generated in the shell, an average electric 
field perpendicular to the magnetic field must exist in the shell rest frame. The electric field is not 
uniform, since the system is charge neutral, but exists only over distances of order the width of a 
filament, with values that are proportional to the magnetic field strength. When the gyroradius is 
large, the electron effectively sees an average electric field as it completes one orbit. This electric 
field gives the electron guiding center a velocity of —U" in the shell rest frame. The effective 
magnetic field strength in the electron rest frame is then 

B" = ^j. (30) 

1 s 

Because the electron travels over many filaments, the orientation of the magnetic field 
changes dramatically over one gyroradius, so that the direction of the electron's velocity vector 
is randomized through a random-walk process rather than through the rotation over one orbit. 
The electron travels the filament width Xf in the time c St". In this time, the angle 9 « St"/r e is 
traveled, where r e is the electron gyroradius. Because the motion is a random walk, the number 
of time intervals St" required to change direction by 2tt is approximately n = 47r 2 /# 2 . The total 
amount of time required to isotropize the motion of the electrons is therefore given by 

„ . „ 4vr 2 r 2 4vr 2 cr 2 , . 

t' so = nAt" = e - = , 31 

CXf LU" Z Xf 

where uj" is the cyclotron frequency. Using equations (19) and (30) to define the magnetic field in 
equation (31), one finds 

C = nAt = • = 11.1 s n is l V J ;= • (32) 



LO, 



p,e,tsm 



f? 3 / 2 ' lsm '-" y/l + 7? 2 + 2nT 



For r" « 1, this timescale is long compared to the two-stream instability for electrons. 



6. Burst Radiation 

Two radiative processes are present within the theory: synchrotron emission and Compton 
scattering. The synchrotron emission will occur isotropically in the electron rest frame describe by 
equation (23). The observed synchrotron radiation will be boosted into the observer's reference 
frame by a factor of T e . Compton scattering of synchrotron radiation by the synchrotron-emitting 
electrons boosts the radiation by another factor of T 2 , because the electrons in this rest frame 
have a characteristic energy of m e c 2 T e . 

The characteristic synchrotron frequency in the shell rest frame is given by 
hu"/m e c 2 = 2B"Fl/3B cr , where B cr = eh/m 2 c 3 . Transforming this into the observer's 
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reference frame and using equations (19) and (30) to remove B", one finds that the characteristic 
synchrotron energy is 

hv s ly/Mhe i i r^r 

:V*nl sm ^r, (33) 



3m yi c 2 V iam (l + Tr 2 + 2T,T)(l + T,T)' 1 ' 

II r 4 

= 2.17 x 1(T 6 v\nf- £ , (35) 

/-3 ^(i + 2 7? r)(l + ?? r) ' v ; 

where terms of order r] and higher were dropped in the last equation. Equation (35) places the 
characteristic energy in the optical band for the given characteristic energies and rjT < 1. As 
rj increases above T , the characteristic observed energy falls as rj 2 , so r] = T defines the 
maximum characteristic frequency for a given value of T. To have a characteristic photon energy 
of m e c 2 at 77 = T^ 1 requires T > 6.93 x 10 4 n^JJ . 



The characteristic energy of the Compton emission after a single scattering is a factor of V 
larger than the synchrotron characteristic energy, so 



hue 2 ^%U (r + q) 5 r m 

: 7 i 2n i Sm -—. , . „ „,2 „ . ^7 > v<30J 



3m 3 e /2 c* ' lsm (1 + V 2 + 2^) 2 (1 + V T) ' 
= 2.17 An^, ^ . (37) 

/ 3 ^( 1 + 2r? r) 2 (i + r/r) v ; 

The characteristic energy of the Compton scattered radiation after a single scattering is at m e c 2 
when r > 1.47 x 10 3 . A second scattering takes the photon in the observer rest frame to 
the GeV energy range. Further scattering does not change the photon energy, because the photon 
energy is of order the characteristic electron energy after the second scattering. 

Each of the radiative components spans a broad range of energies. The low end 
of the synchrotron emission is set by the cyclotron frequency, which is smaller than the 
characteristic synchrotron frequency by a factor of Tg. For r]F = 1, the cyclotron frequency is at 
v ~ 2.68 x 10 s Hz f\i m ^ 2 , so that synchrotron emission extends down to the radio band. An 
important point is that the cyclotron frequency is related to the plasma frequency in the electron 
rest frame as rs y/2ujp e r]Y /2ir (1 + 277F) 1 / 4 , so that they are about equal, and the cyclotron 
photons escape the shell to the observer. The lowest energy of the Compton scattered radiation is 
the cyclotron frequency photons upscattered by Tg , which means that the low end of the Compton 
scattered radiation equals the high end of the synchrotron emission. If most of the energy is in the 
Compton scattered component, and the synchrotron photon number spectrum falls faster than 
v~ 2 , so that most of the energy is released at the low end of the spectrum, then the low end of the 
Compton spectrum will be larger than the high end of the synchrotron spectrum. This implies 
that optical and ultraviolet emission is part of a single smooth continuum that extends through 
the x-ray and gamma-ray bands, and that most of the energy emitted by the burst is released in 
the optical and ultraviolet. 
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The ratio of the synchrotron emission rate to the single-scattering Compton emission 
rate for a single electron is given by P' S y nc /P"i = W'B/W" ynch , where Wg and W" ynch are 
the magnetic field and the synchrotron photon energy densities as measured in the electron 
rest frame. The synchrotron energy density is related to the single electron emission rate by 
Wsy nc = Psy n ch n emisl" /^c, where n" mis is the density of electrons with Lorentz factor T e in the 
electron rest frame. The synchrotron emission rate is given by P's ync = 4aTcTlWg/3, so the ratio 
becomes 

p n , p n Mr/T+l) ( 9m p \^ 

SynJ Cl a T (T + V ff emis {Mnl m ^) ' W 

~ ^ 5Ui3 (r + ?? ) 2 n ^m^27 terms [ R J ■ V») 

The ratio of the synchrotron emission to the rate at which energy carried by the interstellar 
medium flows through the shell is given by 

Knc AmeaTM^nlj^lrj (1 + V T) f erms ( Rq\ 2 . , 



m p u-; Hon ± e 3fm 3 



R 



2 



1.956 xlO- B Xl 7 n| m rJ(r ? r)(l + r ? r)/ emis f§) . (1.1) 



For single scattering Compton cooling, this ratio is 



E'com P _ 2m e alMlnl m TT V {T + r 1 f f emis fR^ ( ^ 



m p c*n ion rl 3 f V R 



= 2.80xl0- 7 ^| 7 4 m r|(r^i^^/ e 2 mis ^) 4 . (43) 

From these equations, one sees that the energy release is very inefficient unless rii sm or V are larger 
than the characteristic values used in the calculations. We discuss this point in §8. 



7. Radiative Timescales 

The cooling time scale for synchrotron emission of a relativistic electron is found by dividing 
the electron energy m e c 2 T e by P'J ync , the emissivity of a single electron. In the electron rest frame, 
the synchrotron cooling timescale is 

f = ? ^ r + 1 ) 2 (44) 

sync Aa T cl 1 n ism T^T + r ] )^lTv ^ +W , 

= 3.76x10^^-^3 2 {f]T / + 1)2 (45) 
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The synchrotron self-Compton cooling timescale is found by multiplying equations (38) and (45) 
together, which gives 



9vr 



9m r 



Comp 



^cfemittriism \ni m M J rjTs (r + r/) 3 y/l + r/ 2 + 2rjT 



R 



1.32 x io 8 S M?C/ir3 3 



far + i^r 



3 r 2 



AY 



(46) 
(47) 



r? (r + V ) 3 ^/1 + V 2 + 2 V T \R J 

These timescales are plotted in Figures 3 and 4 for rij sm = lcm -3 and 10 5 cm~ 3 , with 
(R/Rq) 2 = m p /m e , and nT = 1. One sees that for the higher densities, the Compton cooling 
timescale can fall below the timescale for isotropization. When this occurs, the radiative cooling 
will determine the shape of the electron distribution. The Compton cooling timescale is shorter 
than the isotropization timescale when 



n 



■ism > 1.16 x 10 6 cm 



-3 



V 2T2 (7]T + iy 
(r + r/) 3 



(48) 



8. Observational Consequences 

The theory outlined above has within it two observational selection effects that define lower 
limits on the values of rii sm and T. These selection effects arise because gamma-ray bursts are 
recognized as such through the efficient emission of gamma-rays. 

Because gamma-ray bursts are selected by their gamma-ray emission, one must have a value 
of r\ that is sufficiently small to give gamma-rays through Compton scattering. Because the largest 
value of the photon energy occurs at r/ = 1/T, one can derive a lower limits on T by requiring the 
right hand side of equation (37) be > 1 at this value of rj: 

T 3 > 0.831 n~ s t ■ (49) 

The weak dependence on rij sm implies that for all gamma-ray bursts, the bulk Lorentz factor 
r > 10 3 . Events may occur with smaller F, but these would emit in the optical and ultraviolet 
bands. 

From Figures 6 and 7, one sees that an upper limit on the value of rj is found from equation 
(37) for nF » 1. This limit is 

V < rimax = 1-03 x IO" 3 nf sm rf . (50) 

These upper limits are plotted in Figures 1 and 2. An important aspect of these limits is that the 
timescales associated with the shell thickness are generally < 1 s. For the higher density figure, the 
time scales at R/Rq = 0.1 range from 0.01s to Is, which is consistent with the shortest timescales 
exhibited by gamma-ray bursts. 
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Fig. 6: Characteristic radiation energies. — The characteristic synchrotron emission (dotted lines) 
and synchrotron self-Compton emission (solid lines) energies of photons emitted from the relativistic 
shell as measured by an observer at rest in the interstellar medium for photons emitted along the 
velocity vector of the shell. Values of V = 10 3 (upper curves) and 10 4 (lower curves) are used. The 
density is rjj sm = 1 cm~ 3 . 
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Lower limits on the interstellar medium density are found by requiring that Compton 
scattering efficiently remove energy from the shell. Two conditions must be met: first, Compton 
cooling must dominate synchrotron cooling, and second, the Compton cooling rate must be 
comparable to the rate at which energy is lost as the shell decelerate over the distance R. The 
first of these conditions is derived from equation (39): 

n- > 6 54 AO r f ^ r+1 ) lr3 r -4 ( R V (51) 
n lsm > 0.04 yw 27 j emis ^ r + ?? ^ L s y Ro J ■ K 131 ) 

The second of these conditions is found by equating the right hand side of equation (43) to 
(Ro/R) 3 g, where g < 1 is a measure of efficiency. This gives 

3 

n lsm = 8.22 x 10 4 cm" 3 M^V^hL (^) ' 9 l ■ (52) 

If g is very small in equation (52), then most of the energy lost by the shell to the interstellar 
medium is not radiated away, making the burst dim and unobservable. 

The limits from equations (51) and (52) on rii sm are plotted as functions of T in Figure 8. 
The point to note in this figure is that the density required for high efficiency cooling is very high, 
of order 10 5 cm~ 3 . This is important in explaining why all gamma-ray bursts are have a peak 
value of the vF v curve near 250 keV. The high source density provides a medium for Compton 
attenuation of the burst spectrum. Because the scattering medium is at rest in the galaxy rest 
frame, the value of E p is independent of T, even though the characteristic energy emitted by the 
shell is a strong function of T. As a consequence, if the medium density is low enough to keep 
attenuation from occurring, then the density will be too low to efficiently produce gamma-ray 
emission. In such a circumstance, the shell will loose energy by thermalizing ions and electrons, 
but the thermal energy of the electrons will not be rapidly radiated away, making these particular 
sources invisible. 

An interesting consequence of the density on the theory is that as the density increases, the 
timescale for radiative cooling becomes shorter than the timescale for isotropization of the electron 
distribution. The implications of this is that there should be a coupling of the electron distribution 
to the density of the interstellar medium, with the distribution falling more rapidly, and being 
more anisotropic, at the higher interstellar medium densities. Because the shape of the electron 
distribution determines the shape of the spectrum, one expects the burst spectrum to become 
softer for the higher interstellar medium densities, and therefore for the higher attenuation optical 
depths. This is the case observationally, so this theory may provide an explanation for that one 
characteristic of the Compton attenuation theory. 

The theory as constructed has two natural timescales, one from the thickness of the shell, and 
the second from the deceleration distance R. The shell thickness timescale is given by equation 
(12), and is of order 2.9s for rii sm = 10 5 cm -3 and R/Rq = 1. For the same density and R/Rq = 10, 
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Fig. 8: Limits on density. — The lower limits on rii sm as a function of T. 
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the timescale falls to 0.029s. The deceleration timescale is given by 



t B = JL = 0.203s M$n£fc»' (£) , (53) 



where the definition of Rq in equation (5) has been used. This timescale is less than the shell 
thickness timescale whenever 

— > 8.74 r| 77I3 . (54) 

Because should be of order m p /m e , the two timescales are of the same order. When the two 
timescales are equal, the value is given by 

t R = t shM = 1.77s M^n^T- 5/3 V l 3 . (55) 

For 

n-ism — 10 cm , the timescale is 0.038 s. The timescales associated with both the shell 
and the deceleration distance are sufficiently short to be responsible for the shortest timescales 
observed in gamma-ray bursts. 



9. Summary of Conclusions 

To summarize the theory presented above, a baryonic shell with an ultra-relativistic bulk 
velocity interacts with the interstellar medium through the filamentation and the two-stream 
plasma instabilities. The former instability gives rise to a magnetic field with a strength that is 
far below the equipartition value, and the latter instability heats the electrons to energies that are 
relativistic, but also far below the equipartition value. Neither instability is sufficient to produce 
a shock. Instead, the interstellar medium passes through the shell, so that the region behind the 
shell is not cleared of interstellar medium. The electrons within the shell produce synchrotron 
radiation with a characteristic energy in the observer's rest frame that ranges from the radio to 
the ultraviolet. Synchrotron self-Compton emission by the electrons in the shell produces x-rays 
and gamma-rays in the observer's rest frame, in addition to optical emission. The optical emission 
is dominated by the synchrotron self-Compton component. The timescales associated with the 
shell thickness and the length scale over which the shell decelerates provide a lower limit on the 
burst durations. The burst duration itself would be determined by the complex structure of the 
relativistic wind, since the interstellar medium remains in place, permitting multiple shells to each 
produce gamma-ray emission. 

Two conditions must be met for the interaction between shell and interstellar medium to 
efficiently produce gamma-rays: first, the bulk Lorentz factor must be > 10 3 in order to produce 
radiation above IkeV; second, the number density of the interstellar medium must be greater than 

— 13 

pa 10 6 cm~ 3 r 3 4 , where T 3 is the bulk Lorentz factor in units of 10 3 , for the thermal energy to 
be radiated efficiently. The lower limit on T through the selection effect provides an explanation 
for why the value of T is always sufficiently high to allow the escape of lMeV photons from the 
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gamma-ray burst emission region without the production of an electron-positron plasma from 
photon-photon pair creation and the subsequent thermalization of the radiation. The limit on 
density provides an explanation of why all gamma-ray burst spectra appear to be Compton 
attenuated. Because the limits on T and rii sm are from selection effects in observing the emission 
of gamma-rays, one expects there to be burst events with values of F and rii sm outside these 
limits. For bursts with low V, the bursts radiate at energies below 1 keV. Therefore, one expects a 
population of optical and ultraviolet transients that have no gamma-ray emission. For bursts with 
low density, the radiation of energy is inefficient, so that the bursts are of low intensity. These 
bursts may appear in burst samples through a correlation of burst intensity with the interstellar 
medium density inferred from the Compton attenuation model. 

An aspect of the theory that provides a test is the comparison of instantaneous gamma-ray, 
x-ray, and optical emission to the radio and optical emission. There are two aspects of the 
theory to test. First, one can test whether the broad band spectrum is consistent with being a 
synchrotron spectrum at low frequency and a Compton scattered synchrotron spectrum at higher 
frequency. Second, one can test the consistency of physical parameters in the theory. This last 
is done by comparing the Thomson and photoelectric optical depths found through a fit of the 
Compton attenuation model to the optical attenuation derived under the assumption that the 
unattenuated gamma-ray continuum extends to the optical band. Third, if the optical spectrum 
is sufficient to fix the cyclotron frequency by determining the low energy drop-off of the Compton 
spectrum, one can solve for the value of the Lorenz factor, which will provide a consistency test 
through the lower limit on T. If one can model the x-ray afterglow of a burst as the forward 
scattering of x-rays by dust, then one has additional information about the optical extinction that 
can be used in this comparison. 

The theory of afterglows in this theory has yet to be developed. An important difference from 
the shock theory of afterglows is that the region behind the shell will emit afterglow radiation in 
competition with the radiation from the decelerated shell. Because the evolution of the afterglow 
from the interstellar medium is determined by the broadening of the look-back surface and the 
radiative cooling of the interstellar medium, while the evolution of the shell radiation is determined 
by the decrease in V and the evolution of the thermal structure of the shell, the theory should 
have two distinct afterglow components that produce a complex afterglow behavior. 

Numerical modeling of the plasma processes can lead to additional observational tests of the 
theory. In particular, the numerical modeling of the electron distribution for the regime where the 
electron isotropization timescale exceeds the Compton cooling timescale (Fig. 2) may provide a 
test through the correlation of Thomson optical depth with spectral hardness. One suspects that 
as T increases, the synchrotron spectrum becomes softer, because an electron radiatively cools 
before it isotropizes, making the electron distribution one-dimensional. Because of the lower limit 
on rii sm in inversely related to T, one expects rij sm , and therefore the Thomson optical depth, to 
be smaller for softer spectra. This conjecture requires numerical verification; if it is verified, then 
one can use the correlation of unattenuated spectral index with Thomson optical as a test of the 
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theory. There is already some evidence of such a correlation (Brainerd et al. 1998, Fig. 8a). 

Two theoretical investigations are now warranted. The first is a study of the broad band 
spectrum expected for this theory need to be numerically calculated for a number of electron 
distributions. This study will determine what aspects of the spectrum provide tests of the theory 
without requiring a full understanding of the plasma processes. Such a study can be carried 
out through Monte Carlo simulation, and should include the effects of optical and Compton 
attenuation. The goal is to model the spectrum from the radio to the gamma-ray. The second 
is a study of the plasma processes. This will require the development of plasma codes to study 
the interactions of relativistic beams and the growth of instabilities to the nonlinear regime. Only 
such a study will verify if the analytic conclusions reached above are accurate. Only such a study 
will one determine the value of r\ in terms of the other burst parameters. And only through such 
a study will one obtain model spectra that are dependent just on the bulk Lorentz factor, the 
density of the interstellar medium, and the mass of the relativistic shell. 

The plasma instability theory is capable of explaining the most important features of 
gamma-ray bursts. Further theoretical research is therefore justified, and should lead to strong 
and unambiguous observational tests of the theory. 
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